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O ' ABSTRACT 

<n ; 

(^i • The decay of half of the GRB radio afterglows with long temporal monitoring is 

1} ' significantly slower than at optical frequencies, contrary to what is expected in the 

simplest fireball model. We investigate four ways to decouple the radio and optical 
' decays: an evolving index of the power-law distribution of the shock-accelerated elec- 

trons, the presence of a spectral break between the two domains, a structured outflow, 
£N| • and a long-lived reverse shock contribution to the afterglow emission. For most after- 

' glows, the first scenario cannot accommodate all properties of the afterglow emission. 

CO , If the spectral break of the second scenario is the cooling frequency, then observations 

require that, as the fireball decelerates, the parameter for the minimal electron energy 
decreases and the magnetic field strength is constant or increases. The latter behavior 
seems implausible and is alleviated if the circumburst density increases outward. In 
the framework of the third scenario, the optical afterglow arises from a more energetic, 
narrow outflow core while the radio emission comes from a more extended envelope. 
This scenario is at best marginally consistent with the general properties of the ra- 
Q ' dio and optical afterglow emissions and requires a total electron energy exceeding 

equipartition, thus it does not provide an acceptable solution. In the fourth scenario, 
q ■ it is assumed that the radio afterglow emission arises in the reverse shock propagating 

?-h ' in a steady stream of ejecta, which catch up with the decelerating GRB remnant. This 

c/3 , scenario can accommodate the properties of the afterglows with slow radio decays and 

requires that the injected energy is less than or comparable to the initial fireball energy, 
while other afterglow parameters have reasonable values. We find the reverse-forward 
shock scenario to be the most viable explanation for the shallow decay observed in 
some radio afterglows. For a jet, the transition to a semi-relativistic motion mitigates 
the radio decay, however this scenario would work only when the slower radio decay is 
observed well after the steeper optical fall-off. A structured outflow with a relativistic 
core, yielding a fast decaying optical emission, and a non-relativistic envelope, produc- 
ing a slowly falling-off radio afterglow, is not a viable solution, as it fails to decouple 
the radio and optical decays. 

Key words: gamma-rays: bursts - ISM: jets and outflows - radiation mechanisms: 
non-thermal - shock waves 
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1 INTRODUCTION 

Since their first discovery, the optical afterglows of Gamma-Ray Bursts (GRBs) have displayed a time evolution which is 
consistent with the predictions of the fireball model (Meszaros & Rees 1997). The power-law decay expected for a spherical 
outflow has been first observed in GRB 970228 (Wijers, Rees & Meszaros 1997). A tight collimation of the GRB outflow (i.e. a 
narrow jet) was shown to yield a steepening ("break") of the afterglow light-curve fall-off (Rhoads 1999), which subsequently 
was observed in GRB 990123 (Kulkarni et al. 1999). To date, there are about ten afterglows exhibiting an optical light-curve 
break and a smaller number of sufficiently well monitored afterglows whose light-curves are pure power-laws. In all cases 
where the optical spectral slope was determined, it was found to be consistent (within the fireball/jet model) with the optical 
light-curve decay rate. 

Further tests of the basic afterglow model require observations at earlier times, within the first couple of hours after the 
burst, and at other wavelengths. The number of GRB afterglows with early observations will be greatly enhanced by the 
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Swift satellite. On the latter avenue, the monitoring of about ten X-ray afterglows has shown that their behavior is consistent 
with that observed in the optical, without major surprises for the fireball model. However, it should be noted that X-ray 
observations rarely cover more than one decade in time. Wider time ranges have been reached for about ten radio afterglows, 
half of them exhibiting a decay which is much slower than in the optical, which is a rather puzzling feature for the fireball 
model. That radio light-curves fall-off after about 10 days, strongly indicate that peak of the synchrotron spectrum is redward 
of the radio domain after that time. Then, within the fireball model, one expects nearly equally fast decays in the radio and 
in the optical, contrary to some observations. 

The purpose of this paper is to analyze in more detail whether the fireball/jet model in its simplest form can explain 
the discrepancy between the radio and optical decays observed in some afterglows (§2) and to explore departures from the 
usual assumptions that could accommodate such discrepancies (§3). The expression "simplest fireball model" will be used 
to designate an outflow i) with a uniform angular distribution of the kinetic energy, ii) without energy injection during 
the afterglow phase, Hi) interacting with a circumburst medium whose density changes monotonically with radius, iv) with 
un-evolving and uniform parameters for the electron distribution and magnetic field behind the forward shock. For an easier 
understanding of the following equations, we list in Table 1 the definitions of the quantities which will be encountered more 
often. 



2 RADIO LIGHT-CURVES IN THE SIMPLEST FIREBALL/ JET MODEL 

For a spherical, uniform fireball, or a wide jet whose dynamics is not yet affected by its lateral spreading, the synchrotron 
characteristic frequency corresponding to the lowest energy electrons accelerated by the forward shock (hereafter called 
"injection" frequency) is given by (e.g. Panaitescu & Kumar 2000) 

(^"tel'^lV* CD 

where £ is the fireball kinetic energy per solid angle, td is the observer frame time measured in days, z the burst redshift, 
es is the fractional post-shock energy stored in the magnetic field, and e; parameterizes the lowest energy jim e c 2 of the 
power-law distribution of shock-accelerated electrons, representing the fraction of the energy in the downstream fluid residing 
in electrons if they all had the same random Lorentz factor 7;, i.e. 7;m e c 2 n e = u x rm p c 2 n p , where n e = n p are the electron 
and proton number densities in the post-shock fluid. 

The synchrotron characteristic frequency corresponding to an electron energy for which the radiative losses timescale 
equals the fireball age (called the "cooling" frequency) is 



for a homogeneous circumburst medium of particle density n, and 

, c = 4 x io" ( , + ir 3/2 ( ios^) A - 2 (ifO " /2 {Y + 1)_2 tT Hz ' (3) 



for an r~ 2 wind-like density profile, A* being the ratio of the mass loss rate to speed of the wind, in units of 
lO _8 (M0/yr _1 )(km/s) _1 . In the above equations, Y is the Compton parameter. 

If the afterglow is optically thin at the peak of the spectrum, the synchrotron self-absorption frequency is 

for a homogeneous external medium, and 

-- 3 ^ +I >- 2 "(To4s)" 2 ' S - 4 - ,, (ra)"'(i^)" 5 '" ,S ° H - <» 

for a wind. In the derivation of the above equations it was assumed that v a < Vi < v c . For a high magnetic field parameter 
es, the cooling frequency may be below the injection, in which case v a is given by 

- - » <•+■>-"■ fe) ""' (i^r m"* <>- + •> ^ °- ■ 

for a homogeneous external medium, 

"» = 6 (* + i) 375 (ipfcii) A "' h (ifO 6/5 {Y + 1} ^ 8/5 GHz ' (7) 
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for a wind. At late times, the injection frequency falls below the absorption frequency. In this case, the latter evolves as 
v a oc £~( 3 p+ 2 )/( 2 p+ 8 ) (homogeneous medium) and v a oc f ~ 3 (p+ 2 )/( 2 p+ s ) (wind medium), where p is the exponent of the power- 
law energy distribution of the electrons accelerated at the forward shock: dN/d'y oc 7~ p , 7 being the electron random Lorentz 
factor. 

The flux at the peak frequency (= min{V;, v c }) of the synchrotron spectrum is 

*-»<•+•> ) 2 Ms) & "' m "' ■* • 

for a homogeneous external medium and 

2 ( -^—) A, f-^V) t" 1/2 mJy , (9) 
^ 10 53 ergs J V 10- 3 / d ' w 

for a wind-type medium, where <1l{z) is the luminosity distance. 

For a spreading jet, the time-dependence of the spectral breaks is given by 

v a oc t -1 / 5 for v a < Vi < v c , v a <x. t~ 6 ^ 5 for v a < v c < Vi , Vi oc t~ 2 , z/ c = consf , (10) 

while the synchrotron peak flux evolution is F p oc f _1 . For < i/„, the self-absorption frequency evolution is v a oc 

j-2(p+l)/(p+4)_ 

To assess which radio afterglows exhibit a decay consistent with that seen at optical frequencies, we list in Table 2 the 
expected radio light-curves for spherical and collimated outflows, a homogeneous or wind-like medium, and for the only three 
likely orderings of the break frequencies a few days after the burst. These light-curve indices were derived from the shape of 
the synchrotron spectrum (Sari, Piran & Narayan 1998) and the evolution of the break frequencies and peak flux. Most of 
the cases shown in Table 2 for v a < min{vi, v c } have been presented in other works (e.g. Meszaros , Rees & Wijers 1998; Sari, 
Piran & Narayan 1998; Rhoads 1999; Sari, Piran & Halpern 1999; Chevalier & Li 2000; Panaitescu & Kumar 2000). 

2.1 Comparison with Observations 

As shown in Table 2, the self-absorbed, low frequency emission rises in time in most cases, as the optical thickness of the 
swept-up medium decreases, while the afterglow light-curve at higher frequencies, above both v a and Vi, falls-off. For observing 
frequencies higher than both Vi and v c , the light-curve is £ _ ( 3 p -2 )/ 4 for a fireball (either type of medium) and t~ p for a spreading 
jet. Such behaviors are observed in the radio, optical, and X-ray light-curves of GRB afterglows, as illustrated in Figure 1, 
displaying all afterglows with a good time coverage in the radio. A flat radio light-curve or a rise slower than i 1 ' 4 can be 
noticed for most afterglows until 20-30 days (70 days for 970508), indicating that v a is never well above the radio domain, as 
in this case the radio emission should exhibit a sharper rise. 

The late time radio emissions of afterglows 970508, 000418, and 021004 exhibit the same decay as in the optical, as 
expected when both domains lie above the injection frequency and the source is optically thin in the radio, thus these 
afterglows are readily consistent with the expectations for the simplest fireball model. However, for all other afterglows, the 
radio emission falls-off significantly slower than in the optical. 

If the cooling frequency v c lies above radio (v r ) but below optical (v ), then the difference between the decay slopes of 
the optical (t~ a °) and radio (t _a!r ) light-curves is Aa or — a — a r — 0.25 for a homogeneous medium and Ao? or = —0.25 
for a wind, the former value being close to that observed in 980703. * . This possibility can be tested by noting that v c < Vi 
implies that the optical and X-ray decays have the same slopes. The upper limit at 5 days on the X-ray emission of 980703 
indicates an X-ray slope larger than a x — 1.22 (Figure 1), thus it may be consistent with the optical slope a = 1.37. 

Index differences larger than Aa or = 1/4, such as those seen in the afterglows 991208, 991216, 000301, 000926, and 
010222, require the presence of a stronger spectral break between radio and optical. However, to obtain a falling-off radio 
light-curve, the source must be optically thin at radio frequencies (Table 2), therefore that spectral break may only be i/,-. The 
v a < v r < Vi cases that yield a decaying radio emission are: i) a spreading jet with v a < v r < Vi and v c < Vi, ii) a fireball 
interacting with a wind and v a < v r < v c < Vi, Hi) a spreading jet with v a < v r < Vi < v c . Below, we analyze each case in 
more detail. 



2.2 Spreading jet, v a < v r < Vi < v and v c < Vi 

In this case the radio decay slope is a r = 1, similar to that exhibited by the afterglow 991208 and 000301. Prior to the time 
tb when the jet edge becomes visible (the jet lateral spreading also becomes important at this time if the jet has sharp edges), 
Vi is given by equation (1), while for t > tb, oc t~ 2 . Therefore 

* If electron cooling is dominated by inverse Compton scatterings and electrons are in the slow cooling regime, then v c decreases slower 
for a homogeneous medium and increases faster for a wind than given in equations (2) and (3) (Panaitescu & Kumar 2000), which 
mitigates the fall-off of the optical light-curve and yields a smaller Ao or , thus this case does not accommodate the observations better 



F p = U(z + 1) 3/2 



d L (z 



d L {z = l) 
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implying that, for tb ~ 1 day (as observed), the electron and magnetic field parameters must be close to equipartition, 
ei = es = 0.1, to maintain <; 10 GHz until the last radio measurement. Because v c is constant for a spreading jet, the 
condition is c (100d) < ^i(100d) is equivalent with v c (tb) < Vi(tb)(tb/W0d) . Using equations (2) and (3) with Y 1 (implied 
by e, ~ es), it leads to n il 200 cm" 3 for a homogeneous medium or A* J; 30 for a wind. Then, equations (6) and (7) give 

v a (tb) Z 10 13 t~* /r ' Hz. At t > t b , v a oc i" 6/5 , therefore ^ a (100d) Z 50t~ 2/s GHz, i.e. the source is optically thick at 8 GHz 
at all observing times. Thus the conditions v a < v r and v c < Vi for t Si 100 days are incompatible. 

Furthermore, the condition ^ c (100d) < ^(100d) and equation (11) imply that v c {tb) Si 10 GHz. First, assume that 
Vc(tb) ~ 10 GHz. Then the radio flux F r at tb is the synchrotron peak flux given in equations (8) and (9). The above values 
of e B , n, and A, yield F r (t b ) k, 10 tj^ 2 Jy, decaying as f" 1 at t > t b to F r (100d) Z 100 tVf mjy. For ^a(100d) £ 50 GHz, 
self-absorption reduces the 8 GHz flux to a few mjy, well above the 0.1 mjy radio flux observed at 100 days. The inconsistency 
between the expected and observed radio fluxes is not alleviated if u c {t b ) < 10 GHz (as assumed above) because this requires 
a higher £ or a denser external medium (ee is already maximal), which, due to that F v (tb) oc £ 3 ^ 4 at v c < v < Vi, can only 
increase the radio flux. 

Thus we reach the conclusion that a spreading jet with v a < v r < Vi < v and v c < Vi cannot accommodate the radio and 
optical decays seen in 991208 and 000301. We note that, if the optical emission seen in 991208 after 10 days is not from the 
afterglow but to a supernova contribution (Castro-Tirado et al. 2001), then the radio and optical fall-offs are not measured 
simultaneously, leaving the possibility of a mildly relativistic jet during the radio decay, which yields at v r > Vi a decay slower 
than the t~ p expected for very relativistic jets. However, p S> 1.5 would be required in this case, which is clearly inconsistent 
with the steep t~ 2 ' 5 optical decay of 991208 observed until 10 days. 

2.3 Fireball in wind-like medium, v a < v r < v c < Vi < v a 

In this case a r — 2/3, consistent with the radio decay seen for 991216 and 000926 until 100 days. Although the decrease of Vi 
for a fireball (t~ 3//2 ) is slower than for a spreading jet (t~ 2 ), electron and magnetic field parameters close to equipartition are 
still required in this case to ensure that u c < Vi until 100 days. The reason is that, for a fireball interacting with a wind, v c 
increases as t 1 ^ 2 , thus the Ui/u c ratio decreases equally fast for either a fireball in a wind or for a spreading jet. Consequently, 
this model requires the same dense winds, with A, <; 30, leading to v a = 10 (t/100d)" 8/5 GHz (eq.[7]), i.e. a source that is 
not optically thin at 8 GHz, and to F p ~ 1 (t/WOd)^ 1 '" 2 Jy (eq.[9]), well above the observed radio fluxes. Thus the fireball 
model for the radio decay of 991216 and 000926 encounters the same problems as the jet model for 991208 and 000301. 

2.4 Spreading jet, v a < v r < Vi < v a < v 

The radio decay slope expected for this case, a r = 1/3, is close to that observed in 010222 at 5-100 days. As for the first model 
above, the condition ^(100d) > 10 GHz for a spreading jet with t b ^ 1 day requires that ti and tB are close to equipartition 
(eq.[ll]). 

The other two conditions involved in this hypothesis - v a (tb) < 10 GHz and fi(tb) < v c (tb) ~ set upper bounds on the 
external density n and wind parameter A* . From equations (2)-(5) we find the former condition being more restrictive, leading 
to n ^ 10 cm" 3 for a homogeneous medium and A* Ss 4 t 1 / 2 for a wind. Lower bounds on the same parameters are set by 
that the peak flux at 100 days should exceed that observed at 8 GHz for 010222 (i.e. 0.03 mjy). From equations (8), (9), and 
the F p oc t~ x expected for a spreading jet, we find that n <; 10" 4 tj^ cm" 3 for a homogeneous medium or A* <; 0.02 t~- X J 2 for 
a wind. Given that ^i(100d) <; 10 GHz, n and A* should be close to these lower limits, implying thus a very tenuous external 
medium. 

Substituting these lower limits in equations (2) and (3), it can be shown that v c {tb) ii 5 x 10 16 t 3 / 2 Hz for either type 
of medium, i.e. the cooling frequency is between the optical and X-ray domains. As the cooling frequency is constant for a 
spreading jet, the optical and X-ray decay slopes are expected to be equal. This is consistent with the observations of 010222 
at 0.4-3 days (Figure 1). 

Thus the radio, optical, and X-ray decay slopes of 010222 may be accommodated by a spreading jet with v a < v r < 
Vi < v c < v < v x . We note that numerical calculations of light-curves from spreading jets show that, when the dynamics 
is accurately calculated and the analytical approximation of constant jet radius (Sari, Piran & Halpern 1999) is relaxed, the 
jet light-curve at v < Ui is flat instead of decaying as t" 1 / 3 , thus it may be that, in fact, the radio fall-off of 010222 is not 
consistent with a spreading jet. 



3 BEYOND THE SIMPLEST FIREBALL MODEL 

As argued in the previous section, the radio and optical decay observed for the afterglows 991208, 991216, 000301, 000926, 
and 010222 appear incompatible with the simplest fireball model. For these afterglows, at least one of the assumptions usually 
made in the calculation of analytical light-curves must be relaxed. In scenarios where the radio an optical emissions arise 
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from the same population of electrons, the ratio of the afterglow emission at these frequencies depends on the slope of the 
afterglow spectrum and on the evolution of all spectral break frequencies located between the radio and optical domain. In 
scenarios where the radio and optical emissions arise in different regions of the GRB outflow, the segregation may be either 
angular or radial, i.e. in the directions perpendicular and parallel, respectively, to that of the outflow motion. 

Thus, we consider four ways in which the difference a — a r between the radio and optical decay indices may depart from 
expectations: 

(i) radio and optical arise from the same electron population. 

((a)) there is no spectral break between the radio and the optical domains. Then the variation in time of the slope of the 
synchrotron spectrum, which is determined by that of the electron distribution dN/d-y oc 7~ p , 

((b)) there is a break between radio and optical, whose evolution is different than that expected in the simplest afterglow 
model (§2), 

(ii) radio and optical arise from different electron populations. 

((a)) the anisotropic distribution of the energy per solid angle £ across the outflow structure, so that the radio and optical 
emissions arise from different parts of the outflow, moving at different Lorentz factors (for the same photon-arrival time), 

((b)) the existence of a long-lived flow of slower ejecta which are energized by the reverse shock and dominate the radio 
afterglow emission, while the optical afterglow arises in the forward shock propagating into the circumburst medium. 

We analyze below these possibilities for the five afterglows (991208, 991216, 000301, 000926, and 010222) for which the 
radio and optical decay indices are incompatible with each other in the simplest fireball model. We also discuss the effect of 
non-relativistic motion on the decoupling of the radio and optical light-curve decays. 



3.1 Electron Distribution with Time- Varying Slope 

In a scenario employing an evolving p, the steepening of the optical light-curve decay observed in the afterglows 991216, 
000301, 000926, and 010222 could arise from a non-monotonic evolution of p and not necessarily from a collimated outflow. 
We consider first a spherical outflow and discuss later the case of a jet. 

If the onset of the radio fall-off at t r ~ 10 days is identified with the passage of Vi (which is the most natural explanation) 
then, for a spherical outflow, Vi(t) = v r (t/t r )~ 3 / 2 , independent of the type of external medium. For a power-law density profile 
of the circumburst medium, 

n(r) oc r~ s s < 3 , (12) 

the evolutions of the peak flux F p and cooling frequency are F p oc £~ s /( s ~ 2s ) and v c oc £~( 4 ~ 3s )/( s ~ 2s ) ; respectively. Substituting 
in the afterglow flux at v > v a : 

F - F X J W^r (p - 1)/2 Vi<v<v c , n , 

' 1 ("cM)- (P - 1)/2 (^c)- P/2 Vc<V ' ( 3) 

they lead to a light-curve decay index at frequency v given by 

, , dlnF v 3 , , (\ . v 3 . t \ dp 12 -5s f 4 - 3s 1 ,, 

a ^^-^ = ^ t) + (-2 ln - r + ^rJdih-w^) | + i(4-7)j' (14) 

where the additional last term in the rhs applies only for v c < v. 

Evidently, an evolving index p can yield a pure power-law afterglow light-curve in only one observational domain. Given 
that optical measurements cover the widest time range and usually have the smallest uncertainties, we determine the evolution 
of p by requiring that the optical light-curve is a power-law, as observed, i.e. that the optical decay index given in equation 
(14) does not vary in time. Equation (14) can be recast as 

dp dint , 4 12 -5s f 4 - 3s 1 

where p* = -a + -ry- r < ~ttti r ( \'^ > ) 



p ,- p In^ + fln^ ' -"3 3(4-*) [ 3(4 - s) 

is the electron index that would be required by the observed a if p were constant. Integrating the above equation yields 
In |o + 2 i n z° 

P(t)=P*-(p*-Po) ln ± 2 Z > ( 16 ) 

where p is the value of p at some time t when the optical slope f3 was measured: p = 2/3 + 1 if v a < v c and p — 2f3 if 
Vc < v a . 

There are a few consequences stemming from equations (14), (15), and (16). For p < p*, the index p increases in time, 
while for p > p* it decreases, in either case approaching asymptotically the value p„. The light-curve index at an arbitrary 
frequency v satisfies 

1 , v dp \ , 4 — 3s I ,„„. 



2 v a dlnt \ 4(4 - s) 
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where the sign of the last term in the rhs is "+" for v < v c < v, and "-" for v < v c < v a , otherwise that term is absent. 
Substituting equations (15) and (16) in (17), we obtain 



Therefore, a radio decay significantly shallower than in the optical (a r < a a — 1/4) requires that p < p„, implying an 
increasing p and an optical spectrum that softens in time (df3 /dt > 0). Furthermore, for p < p„, equation (18) gives that 
the afterglow decay in the X-ray should be steeper than in the optical (a x > a ). Also for p < p,, the evolution of ot v , 
dciv/dt oc (p* — p) ln(v /v), shows that the radio decay should steepen in time (da r /dt > 0), while the X-ray decay should 
flatten (da x /dt < 0). 

The above consequences of requiring that a r < a can be used to test the scenario of an evolving index p. Using the 
light-curve characteristics given in Table 3, we find that a radio decay consistent with the observations is obtained only if 
s = and v c < v a , for the GRB afterglows 

(i) 991208, for which equation (18) yields a r (10d) = 0.91 ± 0.14, and av(100d) = 1.17 ± 0.11. The available measurements 
of the optical spectral slope, /?/_ s (3.8ef) = 1.05 ± 0.09 (Castro-Tirado et al. 2001), /3 D (7.5d) = 0.90 ± 0.15 (Djorgovski et al. 
1999), and /3k -«(8.5d) = 0.77 ± 0.14 (Bloom et al. 1998), indicate a hardening of the optical emission that rules out the 
current scenario, in which it is expected that /3 (8d) = 1.10 ± 0.08 

(ii) 991216. For the f3 given in Table 3 for Galactic extinction of E(B — V) = 0.63, the radio decay index is a r (2d) = 0.13± 
0.12, a r (Wd) = 0.52±0.09, and a r (100d) = 0.82±0.07. However, the expected X-ray decay has an index a x (0.1d) = 4.12±0.25 
and a x (2d) = 2.47 ± 0.10, which is too large compared too observations. Similar results are obtained for E(B — V) = 0.46 

(iii) 000301, for which a r (30d) = 0.84 ± 0.08 and a r (100d) = 1.26 ± 0.08. The optical slopes (corrected only for Galactic 
extinction) determined by Rhoads & Fruchter (2001) - /3 K -B(2d) = 0.79 ± 0.03, f3 K - B (3d) = 1.04 ± 0.03, /3 K -s(5d) = 
1.03 ± 0.10, (3 K -R(7.5d) = 0.69 ± 0.10, fas (13d) = 1.44 ± 0.37 - indicate a general softening with A/3 = 0.41 ± 0.38 from 
day 5 to day 13 (i.e. after the optical light-curve break). The expected softening, A/3 = 0.21 ± 0.04, is consistent with the 
observations. 

The evolution of p required to yield the power-law decay of the optical emission of the afterglows 000926 and 010222 leads to 
radio fall-offs that are too steep in comparison with the observations. 

Therefore, the scenario considered here is consistent with the available data in only one case (000301). We note that, a 
collimated outflows leads to a smaller p* , as the effect of collimation is to yield steeper light-curve decays, and thus to larger 
radio decay indices (eq.[18]). Consequently, a jet (spreading laterally or not), is even less able than a spherical outflow in 
explaining the shallow decay observed for the five anomalous afterglows. 

3.2 Time- Varying Electron and Magnetic Field Parameters 

Below the absorption frequency v a , the afterglow light-curve should should rise or be flat, unless v a increased sufficiently 
fast. Such a behavior is not observed, thus the falling-off radio emission rules out that spectral break located between radio 
and optical domains and having a non-standard evolution is v a . If the break were the injection frequency Vi then, to satisfy 
the condition v r < Vi until the last radio measurement requires large values for the es parameter. As shown in §2.2 and 
§2.3, in those cases where the shallow radio decay also requires that v c < Vi, the resulting synchrotron peak flux is too large 
compared to the observed radio flux at about 100 days. As shown in Table 2, the steepest radio fall-off that can be obtained 
for v r < Vi < v c , is t -1 / 3 for a spreading jet. Thus, in the simplest model, the case v r < Vi cannot explain the steep radio 
decays observed for the above five anomalous afterglows, ranging from t~ 2 / 3 to t~ 1,2 . 

The remaining break to be considered is the cooling frequency v c . Its location below the optical domain implies that the 
optical and X-ray light-curves must have the same decay, which is satisfied by the afterglow 000926, but is rather inconsistent 
with 991216 and 010222, for which a x ~a a = 0.21 ± 0.08 and a x -a a = 0.12 ± 0.04, respectively. It is possible that the X-ray 
emission of these two afterglows has a significant contribution from inverse Compton scatterings, which could yield a steeper 
fall-off that at optical frequencies. For 991208 and 000301, there are no X-ray observations to test that v r < v c < Vi. 

The non-standard evolution of v c can be determined from the difference Aa or between the observed radio and optical 
decay indices. Given that the slope /3 of the synchrotron spectrum, F v oc v^ 13 , steepens by 5/3 = 1/2 from v < v c to v > v c , it 
follows that a = a r — (l/2)(dlnv c /d\nt), thus v c (t) oc j- 2 ^-^) The decay indices shown in Figure 1 lead to v c oc t~ is for 
991208, v c oc t~ l a for 991216, v c oc ^ 2 ' 8 for 000301, v c oc ^ 3 ' 2 for 000926, and v c oc t' 17 for 010222. Such steep fall-offs of 
v c exceed by much that allowed in the simplest fireball model, where v c oc i~ x / 2 for a homogeneous medium and v c oc t 1 ' 2 for 
a wind-like medium (eqs.[2] and [3]). For a power-law density profile (eq.[12]) and a spherical outflow, the cooling frequency 
evolution is d In v c /d In t = (3s — 4)/(8 — 2s) . Thus, for n increasing as a small power of the radius (s < 0) , v c decreases roughly 
inversely proportional with t. The fastest fall-off of v c with time is reached for s — » — oo, corresponding to a sharp increase 
in the external density, when v c oc t^ 3 '" 2 , which is slightly slower than that required by the afterglows 991208, 991216, and 
010222, and substantially slower than necessary for 000301 and 000926. 

For this reason, we consider the case where there is a contribution to the evolution of v c arising from a time variation of 
the other afterglow parameters that determine it: es and £ (eqs.[2] and [3]). The magnetic field parameter es could vary in 
time if its value is determined by the Lorentz factor of the fireball, if the shocked external medium mixes with the GRB ejecta 




(18) 
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having a frozen-in magnetic field, or if the circumburst medium is magnetized. The fireball energy £ may decrease in time 
due to radiative losses or may increase if the shock is refreshed through a delayed energy injection caused by a spread in the 
initial Lorentz factor of the ejecta (Rees & Meszaros 1998). Given that we want to explain both the observed radio and optical 
decay indices, a r and a , and that these indices depend on the electron parameter e,, we also allow this parameter to be 
variable. We note that radiative losses couple the variation of £ to those of e, and e B . However, to calculate this dependence, 
a knowledge of the entire electron distribution is necessary. For p > 2, it can be assumed that the distribution extends up 
to arbitrarily high energies, thus its slope and minimal electron energy (parameterized by e^) suffice, however a high energy 
cut-off must exist for p < 2 (as is the case for 991216, 000301 and 010222, if v c < v a ). For this reason, we allow the evolution 
of £ to be independent of that of es and e;. 

The evolution of the afterglow parameters £, ei, and e B can be constrained with the aid of the observed light-curve decay 
indices a r and a . To do this, we chose power-law dependences of the afterglow parameters on the observer time 

£ oc t e , e B oc t b , ei oc t % , (19) 

for ease of calculations and because, in this way, a minimal number of unknowns is introduced. We note that, in a scenario 
where the afterglow parameters are time- varying, a steepening of the afterglow light-curve decay may be due to a non-uniform 
variation of the parameters, thus the optical light-curve breaks observed in the afterglows 991216, 000301, 000926 and 010222 
may not require a tight collimation of the ejecta (i.e. a jet). For this reason, we shall consider both spherical and collimated 
outflows. 



3.2.1 Spherical, Uniform Outflow 

In the assumption of a uniform temperature in the post-shock medium, equal to that set by the shock-jump conditions, the 
fireball dynamics is given by F 2 M oc £, where M oc r 3 ~ s is the mass per solid angle of the energized external medium. 
Together with t oc r/T 2 , it leads to 

log r = const + ^-L-log£-|f^logi, r oc (a) 1/(4 - s) . (20) 

Starting from the basic expressions of the spectral characteristics F p , i>i, and u c (e.g. Sari, Piran & Narayan 1998), equation 
(20) leads to 

log F p = const + |— ^ log £ + i log e B - %^ 2s l0g * ' ^ 

1 1 3 4 35 3 4 3f5 

log Vi = const + - log £ + 2 log a + - log e B - - log t , log v c = const - - - — log £ - - log e B - - - 2g log t . (22) 

For afterglow parameters varying as given in equation (19), equations (13), (21), and (22) yield the following light-curve decay 
indices at radio (z/» < v r < v c ) and optical (y c < v ) frequencies: 

--(^-rhO'+e-w+H 1 *-!"- 1 '-^ < 23 > 

p + 2p-2, 3p-2 
-a = t^-e + t^-b+(p-l)i-^ T - . (24) 

Our choice for v r > Vi is motivated by the uninterrupted power-law decay observed for the five anomalous radio afterglows 
until about 100 days, and the presence of a break in the radio light-curve of the afterglows 991208 (10 days), 000301 (30 days), 
and 000926 (20 days) that, most likely, is caused by the passage of Vi through the radio domain. 

Before using equations (23) and (24) and observations to constrain the evolution of e B and ej, let us determine the 
evolution of £ in the absence of energy injection, i.e. the range of e allowed by radiative losses. The fastest decrease of £ 
thorough radiative losses occurs when the shock-accelerated electrons acquire an energy equal to that in protons and cool 
radiatively faster than adiabatically. In this fully radiative fireball, half of the energy is transferred to the circumburst medium 
during the time a mass equal to that already swept-up is shocked, and is radiated away, i.e. £ oc M^ 1 ■ Taking into account 
that the swept-up mass is Ai oc r 3_s , it yields £ oc r~'- 3 ~ s ' ) . Using r from equation (20), we obtain that, for a fully radiative 
fireball, £ oc T cc t~( 3 ~ s )/( 7 ~ 2s ) (see also Meszaros , Rees & Wijers 1998). Therefore, radiative losses are characterized by 

£ oc t e with - 7 3 < e < . (25) 

For a given value of i, the exponents e and b can be determined from equations (23) and (24) and the electron index 
p = 2(3 required by the slope of the afterglow optical continuum, F v cc v~ fi ° . These slopes, corrected for Galactic extinction 
and possible Small Magellanic Cloud-like dust reddening in the host galaxy (inferred from the curvature of the optical spectrum 
under the assumption that the afterglow intrinsic continuum is a pure power-law), are given in Table 3. In the absence of 
energy injection (i.e. for e given in eq.[25]), we reach the following conclusions for the five anomalous afterglows: 

(i) for s = and e = —3/7 (fully radiative fireball), 6 > 0.9 and i < —1/4. Taking into account that the magnetic field is 
B oc (ne B ) 1/2 r and that T oc i~ 3/7 , it follows that dlogB/dlogt > 0.02 

(ii) for s = and e = (adiabatic fireball), b > 3/4 and i < -2/3. From F cc £~ 3/8 , it results that dlogB/dlogt > 
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(iii) for s = 2 and e = -1/3 (fully radiative fireball), b > 4/3 and i < -1/3. Then V cx t~ 1/3 and r oc t 1/3 (eq.[20]) lead to 
dlogB/dlogt > 

(iv) for s = 2 and e = (adiabatic fireball), b > 1.4 and i < —2/3. From V cc i~ 1//4 and r oc t 1 ^ 2 , it follows that 
dlogB/dlogt > -0.05. 

Therefore, for s = 0, 2 and without energy injection, e» must decrease and es must increase in time. The latter is 
sufficiently strong that leads to a magnetic field which is constant or increasing in time, which is too extreme for plausible 
mechanisms for generating a magnetic field (turbulence and mixing of magnetized ejecta with swept-up medium). For the 
afterglow 991208, the inferred behavior of e» and es imply that, after 10 days, the peak flux F p evolves as t° and t 03 for a 
radiative and an adiabatic outflow, respectively, while the injection frequency z/» evolutions are t~ l s and t~ 2 ' 3 . The resulting 
evolution of F p is rather inconsistent with the radio observations of 991208 at 3-13 days - F p oc £-°- 47 ± - 20 (Galama et al. 
2001), while that of Vi is in agreement with the observed z/j oc t~ 17 7 . 

We note that, without energy injection, a slower increase of es is obtained for an increasing external density. However, 
even in the limit of a sharp increase of the external density (s — ► -co), we obtain that es increases in time. In this case, for 
an adiabatic fireball (e = 0), the index b for the five afterglows satisfies b > 0.1, implying dlogB/dlogt > —0.45, while for 
a radiative outflow (e = —1/2), the resulting b > 0.6 leads to dlogB/dlogt > —0.20, i.e. the magnetic field may decrease 
in time. For the afterglow 991208, the resulting evolution of the peak flux (F p oc t° for a radiative outflow, F p oc t ' 5 for an 
adiabatic outflow) is not consistent with the observations, while that of the injection frequency (yi oc t~ 1,8 and z/» oc t~ 2 ' 7 , 
respectively) favors a radiative fireball. 

If there is energy injection (c > 0), a constant parameter es can be obtained for s = 0. In addition, e» must decrease for 
s = 0, but it must increase for s — 2. If the external medium is homogeneous, the peak flux for the afterglow 991208 should 
increase faster than t 2 ' 2 and the injection frequency should decrease faster than f -5 ' 8 at t ^ 10 days. This is in stark contrast 
with the radio observations of this afterglow (see above). For s = 2, the parameter es increases in time even if there is an 
injection of energy. 



3.2.2 Non-Spreading Jet 

Before investigating a collimated outflow, we consider a structured outflow endowed with a core of a higher energy per solid 
angle £ than outside it. In this case, the tangential motions are slow and the lateral spreading of the core is negligible (Kumar 
& Granot 2003). Because the dynamics of the core remains unchanged, the spectral breaks continue to evolve as given in 
equation (22). However, when the boundary of the core becomes visible to the observer, the increase of the area radiating the 
afterglow emission changes from {r/T) 2 to (rd c ) 2 , where 6 C is the angular opening of the outflow core. This leads to a faster 
decrease of the peak flux corresponding to a multiplicative factor F 2 . From equations (20) and (21), it follows that when the 
core boundary is seen, the peak flux is 

log F p = const + 10 _ ^ log £ + i log e B - g 6 _ ^ log t , (26) 
and the index of the light-curve decay F„ oc t~ a increases by 

5a = 3 " S ~ e (27) 
4 — s 

at any frequency. For an adiabatic fireball without energy injection (e = 0), 5a — 3/4 for a homogeneous medium and 5a = 1/2 
for a wind (Panaitescu, Meszaros & Rees 1998). Energy injection (e > 0), decreases the steepening, while radiative losses 
(e < 0) increase it. For the latter, the effect is small even for a fully radiative afterglow: 5a — 6/7 for s = and 5a = 2/3 for 
s = 2. The optical afterglows 991216, 000301, 000926, and 010222 exhibited a break at around 1 day (Figure 1), the increase 
in the decay index being 5a ~ 0.4, 2.0, 0.7, 0.6 (Table 3), respectively. Therefore the breaks seen in the optical light-curves of 
the afterglows 991216, 000926, and 010222 are consistent with those expected for a structured outflow endowed with a more 
energetic core, while that of the afterglow 000301 is too large and incompatible with a structured outflow. 

As can be seen from equations (23) and (24), an additive term 5a in the rhs of each equation leads to a different solution 
for i but leaves unchanged the value of b. Therefore, in the scenario of a structured outflow without energy injection (i.e. e 
is constrained by eq.[25]), the magnetic field strength must still be constant or increase in time to accommodate the radio 
and optical decay indices of the five anomalous afterglows, for the usual types of external media (s = 0,2). However, in 
contrast to the uniform fireball case, a structured outflow without energy injection leads to evolutions of the peak flux and 
injection frequency that are consistent with the radio observations of the afterglow 991208 (Galama et al. 2001), for either a 
homogeneous or a wind- like circumburst medium. As for a uniform outflow, to obtain a constant eg for a non-spreading jet 
requires energy injection and leads to behaviors of F p and Vi which are in conflict with the radio observations of 991208. 



3.2.3 Spreading Jet 

If the ejecta are confined to a jet with a sharp edge, its lateral spreading is not impeded and the jet dynamics changes from 
that of a spherical outflow when this spreading becomes important. Rhoads (1999) has shown that, for a homogeneous medium 
(s = 0), the radius of a spreading jet increases only logarithmically with time, most of the increase of the swept-up mass 
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being due to the lateral spreading. For an external density increasing with radius (s < 0), the jet radius should increase even 
slower during the sideways expansion phase. However, the approximation of a constant jet radius should become inadequate 
for a decreasing external density (Kumar & Panaitescu 2001). For ease of calculations we make this approximation, which, 
for s < 3, has the effect of removing the dependence of the afterglow flux on the density profile of the external medium. 
The dynamics of a spreading jet is given by (Rhoads 1999) 

r oc E 1/6 t- 1/2 , M oc E 2/3 t , (28) 

where E is the jet energy and M is the mass of swept-up external matter. From here, we obtain the following expressions for 
the spectral characteristics: 

F v oc E^e^t- 1 , Vi oc £ 2/3 e 2 e^V 2 , v c oc E-' 2/3 e B 3/2 t° . (29) 

Substituting the parameter evolution given in equation (19 in the resulting afterglow light-curves, yields the following decay 
indices 

-a r = ?-^-e + ?-^-b+{p-l)i-p, and - a a = ^y^e + ^-^b + (p - l)i - p . (30) 

The dynamics of a fully radiative jet can be derived from equation (28), taking into account that E oc M' 1 (half of 
energy is radiated away while the swept-up mass doubles). Then E oc V oc t~ 3 / 5 . Therefore, in the absence of energy injection, 
radiative losses yield —3/5 < e < 0. From equations (23) and (24), we obtain that the observed decay indices of the five 
anomalous afterglows require 

(i) b > 4/3 for e = -3/5 (fully radiative jet). From B oc e£ 2 T and Y oc t~ 3/5 , it follows that dlogB/dlogt > 0.07 

(ii) b > 1.1 for e = (adiabatic jet). Using T oc t~ 1/2 , it yields dlogB/dlogi > 0.05 

Therefore, for s — 0, 2 and without energy injection, eg must increase in time and, just as for a spherical outflow, this increase 
is sufficiently strong that it leads to a magnetic field that is constant or increases with time. In contrast, the electron parameter 
ti may either decrease or increase. The resulting behavior of the peak flux F p and injection frequency Vi for the afterglow 
991208 are in agreement with the observations if the jet is adiabatic. To obtain a constant es, £j must be decrease faster than 
t~ 2 ' 5 , which, for the afterglow 991208, leads to the same inconsistencies between the evolutions of F p and Vi and observations 
that were mentioned for a spherical outflow (§3.2.1). 



3.3 Structured Outflow 

So far, we have considered two scenarios where the radio and optical afterglows arise from the same population of electrons. 
If there is no break in the synchrotron spectrum between the radio and optical domains and if the slope of the electron energy 
distribution does not change in time, then the discrepancy between the radio and optical decays of the afterglows 991208, 
991216, 000301, and 000926 must be explained by that these emissions arise from different regions of the GRB outflow. The 
key quantity in the calculation of the afterglow emission is the outflow Lorentz factor, which is determined by the energy per 
solid angle and the degree of collimation of each region. Thus, one can imagine that, in a structured outflow, a part of it has 
a higher kinetic energy and a larger Lorentz factor than another part, so that the former radiates more in the optical and less 
in the radio than the latter. 

Within the structured outflow scenario, an optical decay steeper than in the radio can be obtained if the outflow radiating 
in the optical is tightly collimated. Then, for an adiabatic fireball and without energy injection, the optical decay index increases 
by Sa = 3/4 for s — and by Sa — 1/2 for s = 2 at the time h when the boundary of the optically emitting outflow becomes 
visible (eq.[27]). As discussed in §3.2.2, the Sa expected for a structured outflow is consistent with that observed for the 
afterglows 991216, 000926, and 010222, but not for 000301. Furthermore, the near-infrared-optical spectral slope determined 
for the afterglow 000301 by Jensen etal (2001), f} = 0.57±0.02, even the steepest post-break decay index that can be obtained 
with a structured outflow, a — 1.5/3 + 1 = 1.86 ± 0.03, is too shallow compared with that observed. For these reasons, this 
afterglow will not be considered in this section. 

If the radio emitting outflow is quasi-spherical, then the differential decay Aa or — a (t > tt) — a r is equal to Sa if the 
cooling frequency of the jet emission is blueward of the optical and to Sa + 1/4 = 1 (s = 0) or Sa — 1/4 = 1/4 (s = 2) if 
it is below. For the afterglows 991208, 991216, and 010222, Aa or = 0.8 ± 0.1, 0.9 ± 0.1, and 0.84 ± 0.05, respectively, thus 
the differential decay seen between the radio and optical emission of the first three afterglows may be accommodated by an 
optical jet interacting with a homogeneous medium, particularly if the cooling frequency is redward of the optical domain. 
However, the Aa or = 1.6 ± 0.1 of the afterglow 000926 exceeds that allowed in this scenario, thus allowance for a different 
electron index p of the radio and optical outflows must be made. For generality, we will allow for different indices p in the 
radio and optical outflows of all afterglows, thus the only remaining restriction is that of a homogeneous medium for the 
afterglow 000926, where the observed Sa — 0.7 is inconsistent with a wind-like medium. 

Having outlined the general set-up for a scenario that accounts for the observed radio and optical decays (for four of the 
five anomalous afterglows) - a structured outflow containing a core with a higher energy per solid angle, producing the optical 
afterglow, and yielding a light-curve break when its edge becomes visible, and an envelope of lower energy per solid angle, 
from where the radio afterglow arises - we proceed to constrain this scenario by requiring that the optical and radio outflows 
dominate the emission in their respective domain over the entire time interval when the afterglow decay is observed. These 
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constraints are suggested by that, as it softens, the synchrotron emission from the optical core yields a radio contribution 
that falls-off shallower than that from the radio envelope, while the slower deceleration of the radio envelope yields an optical 
emission that decays slower than that from the optical core. 

For simplicity, we consider that the angular distribution of the ejecta energy is uniform in both the core and the envelope 
of the relativistic outflow, i.e. the fireball angular structure is a top-hat function (co-axial uniform outflows). Such a "dual 
fireball" scenario (Frail et al. 2000) yields the maximal decoupling of the radio and optical emissions. 

That the optical light-curve decays from the earliest observation indicates that the line observer-fireball center is within 
the optical jet. In this case, the radio envelope, being less relativistic than the optical core, becomes visible to the observer 
before the edge of the core, thus the radio light-curve break seen in the afterglows 991208 and 000926 at t r = 10 — 20 days, 
well after the time of the optical break, cannot be due to the outflow structure and relativistic effects. Given that for plausible 
afterglow microphysical parameters the injection frequency Ui is expected to cross the radio domain at around 10 days (eq.[l]), 
it is more natural to associate the radio light-curve breaks of 991208 and 000926 (and of all other radio afterglows in Figure 
1) with the passage of through the radio band. 

Assuming that the parameters ti and es are the same everywhere, the larger energy per solid angle within the core 
implies that its injection frequency is higher than for the envelope, therefore crosses the radio domain at a time 
t+ > t r . For tb < t < t+, the radio emission from the optical jet is that for a spherical outflow, F v oc t 1 ^ 2 for s = or F v oc t° 
for s = 2, corrected for the frequency-independent steepening Sa (=3/4 for s = 0, =1/2 for s = 2) produced by seeing the 
jet edge, therefore Fr°^ oc i -1 / 4 for s = and Fr°^ oc t^ 1 ^ 2 for s = 2. These fall-offs are shallower than those observed in 
all radio afterglows (with the exception of 010222 for s = 2), which we attribute to the radio envelope. For t > £+, the decay 
index of the radio emission from the optical jet is equal to that of the optical light-curve or differs by 1/4, thus F r (0) (i > t+) 
decays steeper than the observed radio light-curves. Therefore, the dominance of the radio emission from the radio envelope 
Fr H \t) is ensured at all times if it is brighter than the radio flux from the optical jet at t = t+: Fr R \t + ) > Fr°\t+). From 
t+ > t r , it follows that 



F r <«>(t + ) = F r (t r )(£)~ ar 



(31) 



where F r (t r ) is the observed radio flux at the time t r ~ 10 days when the radio fall-off starts, though any subsequent time 
can be chosen if the onset of the radio decay has been missed. For tb <t <t+, the radio emission from the optical jet falls-off 

as F, 



(O) 



oc t 1 ^ 4 3 ' (see above), therefore 

1/3 



F^(t + ) = F(°\t b ) 



,(o) 



(*6 



(32) 



where F^°\tb) is the peak synchrotron flux from the optical jet at t = t b . This quantity can be determined from the observed 



optical flux F (t b ) which, taking into account that vf\tb) < v Q , is 
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(33) 



where po is the electron index in the optical jet. From that v\°^(t) oc t 3//2 at all times and for any s, the crossing time t+ 



-3/2 



defined by v\ 



(O) 



(£ + ) = v T is t + = tb[vf \tb) / i , r] 2 ^' i ■ Substituting it in equations (31) and (32) and using equation (33), the 

( H )ff ^ TP(°)(+. \ xriolHc. a limit ™i ■ /(°) < 



above condition F r l (t + ) > iv (t+) yields a lower limit on v\ (if,): 
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(34) 



From equation (34) it can be shown that the least stringent constraint on „<°> (tb) is obtained if ui°\t b ) > v a . Physically, 



this can be understood from that, for a given optical flux, the peak flux of the core emission must be larger for vi°\tb) 



than for vi°\tb) 



> v a . A larger peak flux leads to brighter radio emission from the optical jet at t = t+. To compensate, 



a larger t + is required because Fr°\t) decreases in time. Finally, a larger t+ implies a higher vf \tb). As we shall see, the 
higher is the lower limit on „<°> (t b ), the more it becomes a problem for the scenario considered here, thus we restrict now 
our attention to the less constraining case vi°\tb) > v a , for which (3 = (po ~ l)/2. Because the observed slope of the optical 
afterglow is affected by an uncertain extinction in the host galaxy, we choose to determine the index po from the optical 
decay index a (t > tb), which is better constrained by observations than that before the break. When the jet edge is seen, the 
optical light-curve at v a < v c steepens from F D oc t {3p - z),i to F a oc t ip/i for s = 0, and from F a oc ft 3 *- 1 )/ 4 to F D oc t ( - 3p+1)/4 
for s — 2. Thus po = (4/3)a for s = and po = (4q — l)/3 for s = 2. Substituting po in the above expression of £ leads to 
£ = (2/3) (ao — cur) for either type of external medium. 

Table 3 lists the values of the observables that appear in the rhs of equation (34). Substituting them, we obtain for 
v < vf\tb) the following results: 

(i) GRB 991208. For s = 0, po = 2.37 ± 0.12, slope of intrinsic optical spectrum (3o = (po - l)/2 = 0.69 ± 0.06 (requiring 
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A\° at ~ 0.2 to be consistent with the slope determined by Castro- Tirado et al. 2001), vf\t b ) > 8,000 GHz. For s = 2, 
po = 2.04 ± 0.12, (3o = 0.52 ± 0.06 (requiring A^ ost ~ 0.3), uf\t b ) > 240 GHz. 

(ii) GRB 991216. For s = 0, po = 2.12 ± 0.07, f3 = 0.56 ± 0.03 (consistent with the slope measured by Garnavich et al. 
2000 for a Galactic extinction of E(B - V) = 0.63), uf\t b ) > 6,000 GHz. For s = 2, Po = 1.79 ± 0.07, po = 0.40 ± 0.03 
(requiring A^ ost ^ 0.1), vf\t b ) > 340 GHz. 

(iii) GRB 000926. For s = 0, po = 3.07 ± 0.03, /3 = 1.03 ± 0.01 (consistent with the value inferred by Fynbo et al. 2001 
after correction for host extinction of A^ ost ~ 0.2), ^ 0) (t b ) > 5, 000 GHz. 

(iv) GRB 010222. For s = 0, po = 1.69 ± 0.01, Po = 0.35 ± 0.01 (requiring Ay ost ~ 0.2 to be consistent with the slope 
determined by Jha et al. (2001) and Stanek et al. (2001), thus exceeding the host extinction inferred by Galama et al. 2003 
and Lee et al. 2003, Ay° st ^ 0.1), vf\tb) > 730 GHz. For s = 2, Po = 1.36 ± 0.01, po = 0.18 ± 0.01 (requiring A v ost ~ 0.3), 
i>f\t h ) > 30 GHz. 

If the above values of v\°\t b ) are extrapolated back to the time t m i„ of the first optical measurement, the resulting v\°\tmin) 
is redward of the optical domain, consistent with that the optical light-curves decay from the first observation. We note that, 
for a homogeneous medium, the electron index required by the observed radio decay, pr — (4ov + 3)/3, are compatible with 
the above values of po for all afterglows except 000926. However, for a wind-like medium, the index pn — (4a r + 1) /3 is always 
smaller than po- Thus, the latter type of circumburst medium requires different electron indices in the outflow envelope and 
core, which makes this scenario somewhat contrived. 

The above lower limit on ul°\t b ) can be converted into a lower limit on the injection frequency v^ R \t r ) for the radio 
component at the onset of the radio decay. As argued before, we expect that v\ R ^{t T ) = 10 GHz (or Vj R \t r ) < 10 GHz, if 
the beginning of the radio decay occurred before the first observation) , thus a self-consistency test for the structured outflow 
scenario is to require that the lower limit on v\ R ^{t r ) is below the radio domain. From equation (1), we obtain 

S ^(M = 2 l0S ^"2 l0S ^' (35) 

where £^ and £ < -°* > are the energy densities of the radio and optical outflows, respectively. This ratio can be determined from 
the observed radio and optical fluxes (F r (t r ) and F (t b )) and equations (1), (8), and (9). At an observing frequency above Vi 
and below v c , the afterglow emission is F v {t) = F p (v/vi)- f) (eq.[13]), therefore 



Fo(h) F (0) {t 



(O) 



{tb 



(36) 



where p r = (pn — l)/2 is the slope of the radio spectrum, pr being the electron index of the radio outflow, which can be 
determined from the observed decay index: a r = (3pn — 3)/4 for s — and ot r — (Spu — l)/4. By substituting v\ R \t r ) from 
equation (35) and F p from equation (21), we obtain 

. Fr(tr) fS-33 B r \ £™ f S 3p r \ U vf' \t b ) V 

log = \s^Ys + YJ log W) + Ks^Ys + —) log T r +(3rloe + ?° log • (37) 

Extracting the ratio £ ( R > /£(°> from here and substituting it in equation (35), yields 
/8-3£ \ pl H) (t r ) Fr(t r ) /12-5s\ . h v r vf\t b ) 

-. 1" Pr log -77^7 = log j- ' - — log — + P r log — r— h /3 Q log . (38) 

V4-s J v\°\t b ) F (t b ) V 4-s ) b t r v { °\t b ) Vo 

Note from the above equation that v\ R) '(t r ) oc [^ )(t i ,)]( 2 +/ 3 °)/( 2 +^) for s = and v\ R) {t r ) oc [^ 0) (t b )} (1+M/(1+M for 
s = 2. Because v\ R \t r ) increases with increasing v\°\t b ), it follows that the lower limit on ^°\t b ) given in equation (34) 
sets a lower limit on ^ R \t r ). For the afterglow characteristics given in Table 3, and with p r and p o determined from the 
decay indices ov and a , respectively, equation (38) leads to the following results 

(i) GRB 991208: v\ R] \t T ) > 200 GHz for s = and v (R) \t r ) > 9 GHz for s = 2, 

(ii) GRB 991216: vf> \t r ) > 530 GHz for s = and vf> \t r ) > 40 GHz for s = 2, 

(iii) GRB 000926: u\ R] ' (t r ) > 26 GHz for s = 0, 

(iv) GRB 010222: v (R) (t r ) > 130 GHz for s = and vf' \t r ) > 9 GHz for s = 2. 

Therefore, we find that the condition that the radio emission from the envelope outflow is brighter than that from the optical 
jet leads to that v\ R '(t T ) is above the radio domain instead of the v\ (t r ) — v r expected from the decay of the radio 
afterglows after t r . Only for two GRB afterglows (991208 and 010222), the structured outflow scenario may be compatible 
with the observations, provided that s = 2, which, as discussed above, requires different electron indices in the outflow core 
and envelope. 

In the Vc°\tb) < v case, an even higher lower limit on v\ R ^{t T ) is suggested by that v[°\t b ) is now larger. However, 
for a self-consistent treatment, the result given in equation (38) should be re-derived after adjusting equation (36) for the 
vi°\tb) > v case. The additional constraint on vi°\tb) necessary to calculate the upper limit on v\ R ^(t T ) is provided by 
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requiring that the optical emission from the radio envelope does not overshine that from the optical jet at any time when 
the power-law fall-off of the optical afterglow is observed, otherwise a flattening of the optical light-curve is expecte dt . This 
constraint is satisfied at all times if it is fulfilled at the time t max of the latest optical measurement. Similar to the above 
calculation of the radio emission from the optical core and the radio envelope, the condition Fo°\t max ) > Fi R \t max ) leads 
to an upper limit on the cooling frequency ui R \t r ) for the radio envelope emission. From £ W < £(<=>) and Vc x £i/2 for a 
homogeneous medium, it follows that, for this type of external medium, vi R \t r ) is also an upper limit on vi°\t r ), the cooling 
frequency for the optical core emission, which can be easily converted into an upper limit on vf\t b ). 

In the Vc°\tb) < v case, the analytical calculation of vf^it^) is cumbersome and it is best to convert all constraints on 
the afterglow spectral properties into constraints on the afterglow model parameters. The five parameters £^ R \ £ i -°\ n (or A» 
for s = 2), ei, and es for a structured outflow are constrained by three equations (for v\ {t r ) = v r , F r (t r ), and F (t b )) and 
two inequalities (for vi R \t r ) and i>\°^(tb))- Solving for the allowed parameter ranges in all cases above where u\ R \t r ) ~ 10 
GHz, we find that £ {R) < 10 50 ergs, > 1/4, and A, > 20 (for 991208), A* > 150 (for 991216), n £ 10 5 cm" 3 (for 000926), 
and A, > 200 (for 010222). Such large external densities, characteristic for a red supergiant or a molecular cloud, lead to 
a semi-relativistic motion (F < 1.5) and a self-absorbed radio emission at t > t r . More worrisome, the electron parameter 
ei, together with the electron distribution slope pr < 2 required by the observed radio decay index a r , imply total electron 
energies that exceed equipartition with protons. Thus, the marginally viable cases require implausible afterglow parameters. 



3.4 Reverse Shock in Ejecta Outflow 

In this section we consider the possibility that the slowly decaying radio emission observed in some GRB afterglows arises in 
a long-lived reverse shock crossing some incoming, delayed ejecta, while the optical emission is produced in the forward shock 
which energizes the circumburst medium. An injection of ejecta in the fireball, lasting for the entire duration of the radio 
observations (tens to hundreds of days), is unlikely to be due to a very long-lived central engine. Instead we consider that the 
delayed injection is caused by that the ejecta are expelled by the GRB progenitor with a broad range of initial Lorentz factors 
(Rees & Meszaros 1998), the slower ejecta lagging behind the faster part of the outflow and catching- up with it during the 
afterglow phase, as the leading edge of the fireball is decelerated by the circumburst medium. 

Assuming constant parameters and eb, the steepening of the optical light-curve decay must be associated with a tight 
collimation of the ejecta because a sudden reduction of the injected power cannot explain the steepness of the decay after the 
break time tb- For a jet with sharp edges, the slower moving, leading edge of the outflow spreads laterally at t > tb faster than 
the incoming ejecta, which leads to angular non-uniformities in the Lorentz factor of the fireball. To avoid this complication, 
we assume that the jet does not have sharp edges, i.e. it is embedded in an envelope (as considered in §3.3) which prevents 
its lateral spreading, but yielding a negligible contribution to the afterglow emission. 

The decay of the radio emission from the reverse shock is determined by the evolution of and of the peak flux Fp R \ 
which depend on the distribution of energy with Lorentz factor in the ejecta and on the Lorentz factor Ti(t) of the ejecta which 
catch-up with the leading edge of the fireball at some observer time t. For the dynamics of the latter, we consider the simple 
law £ oc t e for the fireball energy. Most likely, the ejecta are expelled over a timescale much shorter than the observer time 
(i£ 0.1 day) when the afterglow observations are made. Because of this, the Lorentz factor Fi of the freely expanding ejecta 
which enter the afterglow fireball at t il 0.1 days is independent on the duration of the ejecta release. From the kinematics of 
the ejecta-fireball catch-up, one obtains that 

r /4 - s\ x / 2 

r=GM) (e<3 - s) - (39) 

where F is the Lorentz factor of the decelerating fireball (eq.[20]). Since the fireball Lorentz factor at 100 days is less than a 
few, the above equation shows that, in order for the reverse shock to last for ~ 100 days, the GRB progenitor must release 
ejecta whose initial Lorentz factor is Fi < 3 -f- 5. We note that the above ratio of Lorentz factors is time-independent. Because 
the Lorentz factor of the reverse shock in the comoving frame of the unshocked ejecta depends on this ratio, it follows that 
the reverse shock Lorentz factor and the typical electron energy behind the reverse shock do not evolve in time. Then, if 
the magnetic parameter es is also the same (which implies that the magnetic field is identical), the reverse shock injection 
frequency evolves as oc vf" 1 /Y 2 , where v[ F ^ is the injection frequency behind the forward shock. 

The equality of es behind both shocks also implies that the ratio of the reverse to forward shock peak fluxes is just the 
ratio of the corresponding number of emitting electrons, Fp^ /Fp F) — Mi/M, where Mi and M are the ejecta mass and 
that of the swept-up circumburst medium, respectively. For simplicity, we assume that the mass of the ejecta energized by 
the reverse shock increases as a power-law in the Lorentz factor of the incoming ejecta: 

Mi(> Ti) oc T- q , q > , (40) 



t Castro-Tirado et al. (2001) note such a behavior in the afterglow 991208, at 20-30 days. Apart from a substantial contribution from 
the radio outflow, fluctuations in the circumburst medium density, an episode of energy injection, or the emergence of an associated 
supernova, are other plausible explanations for the late-time bump in the optical afterglow 991208 
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therefore Mi oc F q (see eq.[39]). Equation (40) is equivalent to the assumption that the energy distribution in the ejecta is 
d£ i /dT l <xr- q . 

Using equations (22) and (26) with £ oc t e , the evolution of the reverse shock injection frequency and peak flux are 

rilog^ (fl) 6-s-e(2-s) dlogF p (fl) 1 3-e-a 
dlogi 8 -2s ' dlogt 2 l e)+ 8 -2s 9 ' 1 ' 

From here it can be shown that the radio flux from the reverse shock, °c Fp (u^ )^ p ~ 1 ^ 2 , evolves as iv oc t~ ar with 

3 — e 6 — s — e{2 — s) , 3 — e — s 

^ = 4(4-,) (p - 1) -^4~7y g ' ^ 

where p is the electron index for the reverse shock, which we assume to be the same as for the forward shock. The exponent e 
which quantifies the amount of injected energy can be estimated from the observed optical light-curve decay index a (t > t b ). 
From equations (22), (13), and (26), it can be shown that the decay index of the forward shock optical light-curve is 

3-e e(16-5s)-s , 3-e e(6 - s) + s - 2 fAn . 
a a = -^- P 4(4 _ s) if vo < v c and a a = — p 2(4 _ g) lf "c < ^ , (43) 

where p is the electron distribution index for the forward shock, which can be determined from the observed optical spectral 
slope (3 . 

Equations (42) and (43) and the observed decay indices (Table 3) allow the determination of the indices e and q. Generally, 
there are three constraints that these indices must satisfy. First, if the injected energy (£;(> r,)) exceeds that existing initially 
(£o), then the exponents q and e are not independent. From £ ~ £i oc T 1 -", £ oc t e and equation (20), one obtains 

(3 " S)g (44) 



7-2s + q 



However, an energy injection that is dynamically unimportant (£i{> Ti) <C £o) does not necessarily imply a negligible injected 
mass: in the e — * limit, e and q are independent quantities. Given that the radio shallow decays are observed to last for 1-2 
decades in time, the constraint given in equation (44) should be used only for e > 1/3, corresponding to at least a two- fold 
increase of the fireball energy. The second constraint is that the light-curve steepening 8a(e) given in equation (27), resulting 
from seeing the core edge, must be larger than that observed in the optical, 8a (Table 3). The reason is that, if measurements 
did not capture the asymptotic early-time optical light-curve index, then the 5a determined through power-law fits to the 
pre- and post-break emission is an underestimation of the true index increase 6a (e). The third constraint is that e must exceed 
the lowest value allowed by radiative losses: e > —(3 — s)/(7 — 3s). 

Imposing the above constraints for the indices e and q, we find that the energy injection scenario is compatible with the 
observations in the following cases: 

(i) GRB 991208, 

• s = and v a < v c (with e = 0.3, q = 3.2), 

• s = and v c < v (e = 0, q — 2.5), 

• s = 2 and v c < v a (e = —0.1, q = 4.0), 

leading to reverse shock peak fluxes evolutions (eq.[41]) which are consistent with those measured by Galama et al. (2003), 
but to injection frequencies whose evolutions are marginally consistent with the observations, 

(ii) GRB 991216, for /3 a = 0.58 

• s = and v < v c (e = 0, q — 3.3), 

• s = and v c < v a (e = 0, q — 6.6), 

while for f3 a = 0.87 

• s — and v < v c (e = 0.3, q = 3.7), 

• s — and v c < v (e = 0, q — 2.9), 

• s = 2 and u c < v a (e = —0.2, q = 4.3), 

(iii) GRB 000926, s = and v a < v c (e = 0, q = 4.2), 

(iv) GRB 010222, 

• s — and v c <v (e — 0.1, q = 3.4), 

• s = 2 and v c < v (e = —0.1, q — 5.2). 

For the slope of the optical continuum determined by Rhoads & Fruchter (2001), the sharp fall-off of the optical emission 
of 000301 is exceeds that allowed by the core edge becoming visible. Alternatively, that steep decay would require the 
outflow energy to decrease faster than allowed by radiative losses. Therefore, the properties of the afterglow 000301 cannot 
be accommodated by the current scenario, and will be left out for the remainder of this section. 

The above values of the e parameter show that the energy of the ejecta reaching the forward shock over a tenfold increase 
in observer time is at most comparable to the initial ejecta energy. This is due to that, for the observed optical spectral 
slope p o , a significant energy injection reduces the optical light-curve steepening expected when the core edge becomes visible 
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(eq. [27]) below that which is observed. For generality, we maintain the parameter e in the following equations, although e ~ 
for the four anomalous afterglows. 

As a further test for the reverse-forward shock scenario, we require that the forward shock emission is dimmer than that 
from the reverse shock emission at radio frequencies and that the latter does not overshine the former in the optical. For this 
test, we make use of equations (22), (26), and (41) for the evolution of the spectral characteristics, noting that the cooling 
frequency is the same for both shocks, as the internal energy density, bulk Lorentz factor, and comoving frame fireball age 
have the same values behind both shocks. 

From equation (13), the reverse shock optical emission at t > t r , when is below the radio domain, is 

FW(t)=F r (i)(^) (P_1)/2 min|l, (^^j ' j where F r (t) = F r (t r ) Q- ) . (45) 

The condition that the slower decaying reverse shock emission does not overshine that from the forward shock, Fg F \t) = 
F (tb)(t/~tb)~ a ° , until the last optical measurement (t max ), leads to an upper limit on the cooling frequency: 

log < 2 log + (p _ 1} i og V JL + i_|£ ( i + e) , og *22i + 2 a r log ^ + 2a Q log - . (46) 

As for the case of a structured outflow (§3.3), the radio emission F^ from the reverse shock is not overshined by that 
from the forward shock, Fr F \ \t Fr R \t+) > F^ (t+), where t+ is the time when v\ F ^ crosses the radio domain. This condition 
leads to the same lower limit for the forward shock injection frequency (y\ F \t b )) as in equation (34), except that now 



C(e)^ + 



2 3-e 



3-e(13-4s) 
3(4 - s) 



\ > < 47 > 



and the cooling frequency is that given in equation (46). The lower limit on v[ F \tb) implies a lower limit on the forward 
shock injection frequency at the time t m i n of the first optical observation, v[ F \tmin) = v i i F \t b )(t b /t ra i n Y i ^ £ ^ 2 , which must 
be redward of the optical domain, to explain the decay of the optical light-curve. This requirement provides a non-trivial test 
for the reverse-forward shock scenario, as there are cases where, even for the highest injection frequency allowed by equation 
(1), the peak flux of the forward shock emission at the when vf^ reaches the radio domain is brighter than the observed radio 
flux (attributed to the reverse shock). The above cases that satisfy this constraint are: 

(i) GRB 991208 

• s = 0, v (F) (t b ) > 3 x 10 12 Hz, v c (t b ) < v , 

• s = 2, u\ F) {t b ) > 3 x 10 11 Hz, v c (t b ) < 3 x 10 14 Hz, 

(ii) GRB 991216, for f3 = 0.58 

• s = 0, vf\tb) > 10 13 Hz, v c {t b ) > is , 

• s = 2, v^itb) > 6 x 10 12 Hz, u c (t b ) > v , 

while for f3 a = 0.87, s = 0, v {F) \t h ) > 2 x 10 12 Hz, v c (t b ) < v , 

(iii) GRB 000926, s = 0, v\ F) {t b ) > 4 x 10 12 Hz, u c {t b ) > v a , 

(iv) GRB 010222, 

• s = 0, f\ F) {t b ) > 7 x 10 12 Hz, u c {t b ) < 6 x 10 13 Hz, 

• s = 2, v\ F) {t b ) > 10 13 Hz, v c {t b ) < 50 GHz. 

Finally, the reverse-forward shock scenario must satisfy two other constraints: the observed optical flux at the break time, 
F (t b ), and the passage through the radio domain of the reverse shock frequency at time t r when the radio light-curve 
begins to fall-off. The former determines the peak flux of the forward shock emission: 



(p-l)/2 ( / \ 1/2 



F^>(t b )=F (t b ) 



Vo 



max ■ 



Vc{tb) 



(48) 



The latter constrains the fireball Lorentz factor, as following. If the microphysical parameters e; and es have the same values 
behind both shocks, then the ratio of the reverse shock to the forward shock injection frequencies is 

fr'-J-V v, tv i, (r*-r) 2 5 + e- s fAa . 

is the Lorentz factor of the shocked ejecta in the frame of the incoming, undecelerated ejecta. Note that for e = 0, V = 1.25 
for s — and T' = 1.06 for s = 2, therefore the reverse shock crossing the ejecta is semi-relativistic. Taking into account the 
evolution of the reverse shock injection frequency (eq.[41]) and equation (49), the condition vl R \t r ) = v r is equivalent to 

r(t„) - ( - 5 + e -^ ('A) 16 - 4S (50) 

l[tb) - 4(4- a )(e+l) \tj v r ■ (50) 
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With the aid of equations (34), (46), (47), (48), and (50), one can determine u^ F \tt), f c (tb), Fp F \tt), and r(£;,) from 
afterglow observations, which can be turned into four constraints for the four basic afterglow parameters £, n (or A* for a 
wind), e B , and ej, after using equations (1), (2), (3), (8), (9), and the expressions for the fireball Lorentz factor (e.g. Panaitescu 
& Kumar 2000): 

rW = 8(^i) 3/8 ( I ^) 1/8 ( T ^z3)" 1/8 ^ 3/8 for s = o, r W = i5(, + i)V 4 (_|_) 1/4 j4 - 1 /4 t - 1 /4 for s = 2 . (51) 

From here, the allowed ranges of each parameter can be used as a test of the reverse-forward shock scenario by requiring 
that e» and e B are below unity, and that £ and n (or A,) have plausible values. For the cases outlined above, involving a 
homogeneous medium, we find the following constraints on the afterglow parameters 

(i) GRB 991208, £ < 10 51 ergs/sr, n € (0.01, 10) cm" 3 , u e (0.05,0.2), e B > 0.04, 

(ii) GRB 991216, £ < 2 x 10 52 ergs/sr, n < 30 cm" 3 , e; £ (0.05,0.2), e B > 0.02, 
(hi) GRB 000926, £ < 6 x 10 51 ergs/sr, n < 50cm~ 3 , e; € (0.06,0.3), e B > 0.05, 

(iv) GRB 010222, £ < 3 x 10 50 ergs/sr, n G (0.05, 3) cm" 3 , u £ (0.03,0.07), e B > 0.1. 

Compared to the afterglow parameters found by us (Panaitescu & Kumar 2002) through fitting the radio, optical and 
X-ray light-curves of these afterglows with the forward shock model (i.e. without any contribution from the reverse shock), 
the above kinetic energies per solid angle for the reverse-forward shock model are slightly lower, the magnetic field parameter 
is somewhat larger, while the circumburst densities and electron energy parameter are similar. The largest discrepancy occurs 
for the afterglow 010222, for which radio data were not available at the time when our afterglow fits were done. 

For a wind- like medium, the required kinetic energies are below 10 50 ergs/sr, the microphysical parameters are slightly 
larger, and the wind density corresponds to A* £ (0.1, 1). We note that kinetic energies under 10 ergs/sr require an extremely 
efficient GRB mechanism to yield the observed 7-ray outputs, ranging from 10 52 to 5 x 10 52 ergs/sr (Bloom, Frail & Sari 
2001) for the above four afterglows, thus a wind-like medium is somewhat less favoured. 



3.5 Non-Relativistic Outflows 

The scenarios considered so far involved relativistic fireballs. For either a spherical outflow or a jet with sharp edges, the 
afterglow light-curves asymptotically reach the same decay as the outflow becomes non-relativistic because the lateral spreading 
of the jet diminishes and the entire emitting surface is visible in both cases. As shown by Vietri & Waxman (2000), for a jet, 
the transition to the non-relativistic regime mitigates the light-curve decay, while for a fireball it leads to a steeper decay. 
This may suggest that a slow radio decay following a fast optical fall-off may be explained by a jet becoming non-relativistic 
(Frail et al. 2004). 

In the above scenario, the optical light-curve is F oc t~ p (relativistic jet). From the basic equations for dynamics and 
evolution of spectral parameters, the radio flux from a non-relativistic fireball is found to be 

F v%<v<Vc oc t - 3 (5p- 7 )/«> for s = ; F„ ( <„<„ e cc r (7p - s)/6 for s = 2 . (52) 

For the electron distribution index implied by the optical light-curve index a (Table 2) , the resulting radio decay is significantly 
steeper than observed for the afterglows 000301 and 000926, for either type of medium. For the afterglows 991208, 991216, 
and 010222, a homogeneous medium leads to radio decays shallower than observed, while a wind medium yields radio fall-offs 
that are too steep. Since the non-relativistic transition flattens the decay of the emission from a jet, it can be argued that a 
semi-relativistic motion leads to a milder flattening than expected for the non-relativistic regime, becoming thus compatible 
with the afterglows 991208, 991216, and 010222, for which the expected asymptotic radio decay is too shallow in the case of 
a homogeneous medium. 

A severe problem with the scenario of a relativistic optical and semi-relativistic radio afterglow is that it requires the 
shallow radio decay to be seen after the steeper optical fall-off. Therefore this scenario cannot be at work in the afterglows 
991216 and 010222, whose radio and optical observations overlap in time over one and almost two decades, respectively. For 
the afterglows 991208, 000301, and 000926, the overlap is short-lived, still it is rather unlikely that, during the short time 
interval from the last optical measurements and untill the onset of the radio decay, a relativistic jet decelerating as F oc t^ 1 ^ 2 
can slow down to a semi-relativistic motion. 

We note that a scenario where the fireball becomes non-relativistic before the onset of the optical fall-off cannot explain 
the observed index difference Aa or . Above the cooling frequency, the light-curve for a non-relativistic afterglow is 

F Vc<v oc t- 3(p - 2)/4 for s = ; F„ c< „ oc r (7p - 8)/e for s = 2 . (53) 

Using equation (52), it follows that, when the cooling frequency is between the radio and optical domains, Aa or — 0.1 for a 
homogeneous medium and Aa or = —0.5 for a wind, neither of which is compatible with the observations of the five anomalous 
afterglows. 

As the non-relativistic motion of a uniform GRB remnant seems unable to accommodate the five anomalous afterglows, 
we turn our attention to a structured outflow. We consider again the dual outflow scenario, which allows a maximal decoupling 
of the radio emission, arising from the non-relativistic, uniform outflow envelope, from the optical emission produced by a 
relativistic, uniform core. Equation (51) shows that, for the envelope to become non-relativistic at t nr ~ 10 days after the 
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burst, the external density must ben~ 10 (£r/W ) cm or the GRB progenitor must have a wind with A* ~ 10 4 (£V10 53 ) 
(i.e. much larger than for any known type of massive star), where £r is the kinetic energy per solid angle in the radio envelope. 

Substituting in equations (1), (4), and (5), such high densities lead to an optically thick afterglow at the peak of the 
synchrotron spectrum at t = t nr , therefore the self-absorption frequency is given by v a = Virf^ p+4 \ instead of equations (4)- 
(5), where r» = (5e£) / (Bjf) is the optical thickness to self-absorption at the injection frequency, Sanr being the electron 
column density of the swept-up circumburst medium, and e the electron charge. For p ~ 2 and a homogeneous medium with 
the high density required by the early onset of the non-relativistic regime, the resulting self-absorption frequency is 

If the radio domain were self-absorbed at t > t nr , the radio emission would rise as F r oc t 1 ' 1 (for p = 2; Frail, Waxman 
& Kulkarni 2000) during the non-relativistic phase. Therefore the observed radio fall-off implies that the radio emission from 
the envelope is optically thin at t = t nr . Then, equation (54) leads to £ R < 2 x 10 50 (6i/0.03)~ 1/2 (e s /10~ 3 )~ 1/2 erg/sr. 
For such a low energy per solid angle in the envelope, the condition for non-relativistic motion at t = 10 days becomes 
n ~ 300 (6i/0.03)~ 1 / 2 (eB/ / 10~ 3 )~ 1,/2 cm -3 . Making use of equation (51), for this circumburst density, the optical core must 
have an energy per solid angle £o ii 10 54 (ei/0.03) 1//2 (cb/10 -3 ) 1 ' 2 erg/sr for the core to remain relativistic untill at least 10 
days and accommodate the observed steep optical decay. Assuming the same microphysical parameters in the entire outflow, 
that F p oc £ (eq. [8]) implies that the peak flux of the core emission is much larger than that of the envelope emission. As 
it softens, the core emission should become brighter at radio frequencies than the envelope, therefore this scenario cannot 
decouple the radio and optical emissions and cannot account for the five anomalous afterglows. 



4 CONCLUSIONS 

As shown in Figure 1, the power-law decays of the radio and optical emissions of the GRB afterglows 970508, 980703, 000418, 
and 021004 are similar, as is expected in the simplest form of the relativistic fireball model (adiabatic dynamics, uniform 
outflow, no energy injection, constant microphysical parameters). The shallow decay of the radio emission of the afterglow 
010222 is only marginally compatible with that from a spreading jet but requires some extreme parameters (§2.4). However, 
the shallow radio decay and the steep optical fall-off of the afterglows 991208, 991216, 000301, and 000926, cannot be explained 
within the simplest afterglow model (see also Frail et al. 2004). The reason is that, for this model to explain the observations, 
the injection frequency i/j should be between the radio and optical domains until about 100 days, which, as shown in §2.2 and 
§2.3, leads to a radio emission much brighter than observed. 

We are thus forced to consider scenarios that decouple the radio and optical decay indices. Assuming that the radio 
and optical emissions arise from the same population of electrons energized by the forward shock, we have investigated two 
scenarios involving microphysical parameters that vary in time. Alternatively, the low and high frequency emissions may arise 
in different parts of the fireball if the outflow has an angular structure or if there is a long-lived reverse shock. 

In the former category, we considered the case where there is no spectral break between the radio and optical domains, 
requiring a variable slope of the synchrotron spectrum, i.e. a time- varying slope of the electron power-law energy distribution 
dN/d-y oc 7~ p . The evolution of p can be determined (§3.1) by requiring that the resulting optical light-curve is a power-law 
in time, as the longest time coverage with the least measurement errors is achieved in this domain. For the five afterglows 
with shallow radio decays, we obtain that p should increase in time, which implies that the optical spectrum should soften 
in time, the radio light-curves decay should steepen (d 2 F r /dt 2 < 0), while X-ray light-curves should flatten (d 2 F x /dt 2 > 0). 
With the exception of the afterglow 000301, the resulting radio and X-ray light-curve decays and the optical spectral slope 
are in conflict with the observations. 

Also in the first category, we investigated a scenario where the cooling frequency v c is between the radio and optical 
domains. In this case the optical and X-ray decays should to be identical, as for the afterglows 000926 and 010222, but 
inconsistent with the observations for 991216. To explain the shallowness of the radio decays and the steepness of the optical 
light-curves, v c must decrease faster than usually expected (y c oc t^ 1 ^ 2 for a homogeneous medium), which requires an 
increasing magnetic field parameter es. The evolution of es and that of the parameter for minimal electron energy were 
determined from the observed radio and optical decay indices (§3.2). We found that, if there is no energy injection, for any 

geometry of the ejecta (spherical or collimated) and any radiative regime, ti must decrease and es must increase with timet. 
For both a homogeneous or a wind-like medium, the increase of es is so strong that the magnetic field B is either constant or 
increases in time, which is not consistent with a field arising in a turbulent, decelerating flow, but could be compatible with a 
GRB occurring in a pulsar bubble (Konigl & Granot 2002), where the fireball interacts with a magnetized medium, provided 
that a supernova occurred years before the burst. If there is energy injection, solutions with constant es are possible; however, 
this scenario leads to evolutions of the peak flux F p and injection frequency Vi which are in contradiction with the behaviors 
inferred by Galama et al. (2001) for the afterglow 991208 from its radio observations. A magnetic field strength that decreases 

t By modeling the broadband emission of four GRB afterglows with evolving eg, Yost et al. (2003) have found that good fits can be 
obtained for es decreasing slower than t~ 0A or increasing slower than t 3 / 4 
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in time can obtained if the density of the circumburst medium increases with radius. In this case, the resulting evolutions of 
F p and vi are more consistent with the observations of 991208 if the outflow is collimated. 

We did not considered a scenario where all microphysical parameters (ej, es, p) evolve in time, such that the difference 
Q — ct r observed between the radio and optical decay indices is due both to a steepening electron distribution and a cooling 
frequency decreasing faster than usually. The increase of the electron index p would reduce the required increase of es, though 
a significant evolution of the former would probably lead to X-ray light-curve decays that are too steep compared with those 
observed for the afterglows 000926 and 010222. 

In the second type of scenarios, we have considered a dual outflow endowed with a uniform core of higher energy per 
solid angle, emitting the optical afterglow, surrounded by a uniform envelope of a lower energy density, releasing the radio 
afterglow. This scenario allows a better decoupling of the radio and optical emissions than for an outflow with a smooth 
angular structure, and finds support in the steepenings ("breaks") observed in the decay of the optical afterglows 991216, 
000301, 000926, and 010222, at around 1 day after the burst. Such breaks are expected when the boundary of the optical core 
becomes visible to the observer. Another possibility is that the low frequency afterglow emission arises in the semi-relativistic 
reverse shock occurring when there is a continuous inflow of ejecta into the fireball, while the higher frequency emission is (as 
usually) attributed to the forward shock. To explain the observed steepening of the optical decay, we have maintained the 
premise of a structured outflow, but we have ignored the radiation emitted by the envelope. 

Within the framework of these two scenarios, we required that the softening of the emission from the optically radiating 
region does not make its radio contribution overshine the low frequency emission from the radio emitting region. This require- 
ment sets a lower bound (eq.[34]) on the injection frequency for the optical component at the time U (~ 1 day) when 
the optical light-curve break is observed. At the same time, we required that the shallower decaying emission from the radio 
component does not yield an optical contribution exceeding that arising from the optical region. This condition leads to an 
upper limit on the cooling frequency v c (tb) for the radio outflow (eq.[46]). 

For the anisotropic outflow scenario (§3.3), a stringent test is provided by the lower limit (eq.[38]) on the injection 
frequency of the radio emitting envelope which can be derived from that of the optical core. Since the onset of the radio 
decay is due to the passage of the former injection frequency through the radio domain, the lower limit on it should be below 
10 GHz. For the four afterglows (991208, 991216, 000926, 010222) with slowly decaying radio emission and optical breaks 
compatible with a structured outflow, this condition is at best marginally satisfied. In those marginal cases, total electron 
energies exceeding equipartition with protons and very dense external media, leading to a self-absorbed radio emission and 
semi-relativistic fireball. Therefore, the structured outflow scenario cannot account for all the radio and optical properties of 
the above four afterglows. 

In the reverse-forward shock scenario (§3.4), we found that the sets of parameters describing the energy and mass injection 
which are self-consistent and which can explain the magnitude of the break observed in optical light-curves, correspond to a 
modest increase in the fireball energy, by less than a factor 2. To explain the shallowness of the radio decay, the incoming 
ejecta must have an energy distribution which increases rapidly with decreasing Lorentz factor: dEijdVi oc F { ( 2 o • 4 ' for a 
homogeneous medium and oc r i ' 4 ' for a wind-like medium. An important constraint is provided by that the afterglow 
optical light-curves decay from the time t m i n of the first measurement, which requires that the lower bound on the forward 
shock injection frequency is redward of the optical domain at tmin- This lower limit, the upper limit on the cooling frequency 
(see above), the observed optical flux, and the condition that the reverse shock injection frequency crosses the radio domain 
when the radio decay starts, provide four constraints for the fireball kinetic energy per solid angle (£), external particle 
density (n), and microphysical parameters ej, and e_g- We found that, to explain the properties of the radio and optical 
light-curves of the afterglows 991208, 991216, 000926 and 010222, the reverse-forward shock scenario requires total electron 
energy and magnetic field parameters which are close to equipartition, circumburst particle densities between 10~ 3 cm -3 and 
10 cm~ 3 , and low fireball kinetic energies (upper limits ranging from 3 x 10 50 to 3 x 10 52 ergs/sr), which require a very efficient 
GRB mechanism. These parameters are roughly consistent with those inferred through modeling of the afterglow broadband 
emission with the external shock only (Panaitescu & Kumar 2002). Overall, the reverse-forward shock scenario is the most 
promising of all scenarios discussed here in accommodating the shallow radio decay observed for the five anomalous GRB 
afterglows. 

A possibility not discussed so far is that the anomalous radio/optical afterglows are due to a new type of spectral break 
located between the radio and optical domains. Such a feature has been used to explain the broadband emission of the 
afterglows 991208 and 000301 (Li & Chevalier 2001, Panaitescu 2001). The passage of this break frequency (z/*) through the 
optical domain yields a softening of the optical spectral slope simultaneous to the light-curve break, as seems to have been 
the case for the afterglow 000301. To determine the evolution of i/« requires a measurement of the optical spectral slope v„ is 
redward of the optical domain. Such measurements exist only for the afterglows 991208 and 991216, for which we find that 
v* tx t~ 2 for a homogeneous medium and v, cx t^ 1 for a wind-like medium, assuming that the ejecta are spherical. These 
scalings correspond to a break in the electron distribution at energy 7*m e c 2 satisfying 7* oc r 3 ' 5 and 7, oc T 2 ^ 3 , respectively, 
where F is the fireball Lorentz factor. While the former scaling has no obvious interpretation, the latter is not far from the 
7* oc T expected for p < 2, if the total electron energy is assumed to be a constant fraction of the post-shock energy. 
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TABLE 1. Definitions of quantities pertaining to afterglow dynamics and radiation 



Notation Definition 



v a synchrotron self-absorption frequency 

Vi injection frequency (synchrotron characteristic frequency for the typical electron energy) 

v c cooling frequency (synchrotron characteristic frequency for electrons whose radiative cooling timescale equals the fireball age) 

F p flux at the peak of the synchrotron spectrum F v 

a index of the afterglow power-law light-curve F„ oc t~ a 

Sa increase a (t > tb) — a (t < tb) of the optical decay index at tb 

(3 slope of the afterglow power-law spectrum F v oc v" 13 

£ fireball kinetic energy per solid angle, '°- ) for the optical core, ^ for the radio envelope 

p index of power-law electron energy distribution, dNe/d'y oc 7~ p 

€i parameter for the minimum electron energy. For p > 2, the electron total energy is a fraction jE^ti of the post-shock energy 

es fraction of post-shock energy imparted to the magnetic field 

s index of circumburst medium particle density, n(r) oc r~ s 

f r , v typical frequency for the radio and optical domains, respectively 

superscript for the radio emitting outflow envelope or reverse shock 

superscript for the optically emitting outflow core 

superscript for the forward shock 

t m in,tmax epoch of first and last optical measurement, respectively 

tt time when the optical light-curve decay steepens 

F (tb) optical flux at the break-time 

t r time when the radio light-curve begins to fall-off 

F r (t r ) radio flux at onset of radio decay 
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TABLE 2. Temporal scaling of the synchrotron flux from an adiabatic fireball at an 

OBSERVING FREQUENCY V FOR THE THREE ORDERINGS OF THE SYNCHROTRON BREAK FREQUENCIES V a , 
AND V c THAT MAY OCCUR AT 1-100 DAYS AFTER THE BURST 







fireball, 


n oc r 




fireball, 


n oc 


r -2 




jet 




case: 


(1) 


(2) 


(3) 


(1) 


(2) 
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(1) 
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Va <V C < Vi 


t 1 


t l/6 


t -l/4 
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t- 2 /3 
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Vi < V a < V c 


t l/2 


t 5/4 


t -(3p-3)/4 
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•(3p-l)/4 
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(1) : below the lowest spectral break 

(2) : v above the lowest break but below the mid break, 

(3) : v above the mid spectral break but below the highest 



TABLE 3. Parameters of radio, optical, and X-ray light-curves (see Table 1 for defi- 
nitions). Uncertainties of decay indices a are given in Figure 1. Times in days, fluxes in 
mjy. 



GRB a r a f a x t a (i % t b F (t b ) 8a a t r F r (t r ) 

991208 0.97 1.78 4 1.05±0.09 a <2 > 0.1 10 1 

991216 0.69 1.59 1.82 2 0.5$/Q.S7 b 1 0.2 0.37 ± 0.06 <2 >1 

000301 1.26 2.68 3 0.57 ± 0.02 c 5 0.015 2.03 ± 0.08 30 0.3 

000926 0.68 2.30 2.19 1 1.00±0.18 d 1.5 0.03 0.71 ±0.06 20 0.4 

010222 0.43 1.27 1.39 0.2 0.78 ± 0.02 e 0.6 0.06 0.60 ± 0.06 1 0.2 



f optical light-curve decay index measured after the break, i.e. at t > t b 
| slope of optical spectrum determined at epoch t 
a from Castro-Tirado ct al. (2001) 

b 0.58 ± 0.08 for Galactic reddening of E(B - V) = 0.63, 0.87 ± 0.08 for E(B - V) = 0.46 
(Garnavich ct al. 2000, Halpern et al. 2000) 

c after correction for host extinction of Ay = 0.14 ± 0.01 (Jensen et al. 2001) 
d after correction for host extinction of Ay = 0.18 ± 0.06 (Fynbo et al. 2001) 
e after correction for host extinction of Ay ~ 0.1 (Lee et al. 2001, Galama et al. 2003) 
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Figure 1. Radio (8 GHz), Optical (5 X 10 14 Hz), and X-ray (5 keV) light -curves of all GRB afterglows with well-monitored light-curves 
in the two lower frequency domains. Power-law fits are indicated for each frequency range, errors are given in parentheses. If the last two 
optical measurements for 991208 are excluded, then the optical decay is £- 2 - 4 0±0.l6 Some X-ray light-curves have been shifted upward 
by the indicated factors. Note that power-law decay slope of the radio emission of 991208, 991216, 000301, 000926, and 010222 differs 
from that of the corresponding optical fall-off slope by more than 1/4, contrary to what is expected in the simplest fireball model at late 
times, when the injection frequency should be redward of the radio domain. 
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